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ABSTRACT 

We present the hrst resnlts of a deep nnbiased submillimetre survey carried 
out at 450 and 850/xm. We detected 12 sources at 850 /im, giving a surface 
density of sources with Sgso/^m > 2.8 mJy of 0.49 ± 0.16 arcmin”^. The sources 
constitute 20-30% of the background radiation at 850/rm and thus a signihcant 
fraction of the entire background radiation produced by stars. This implies, 
through the connection between metallicity and background radiation, that a 
signihcant fraction of all the stars that have ever been formed were formed in 
objects like those detected here. The combination of their large contribution to 
the background radiation and their extreme bolometric luminosities make these 
objects excellent candidates for being proto-ellipticals. Optical astronomers 
have recently shown that the t/R-luminosity density of the universe increases 
by a factor of ~10 between z = 0 and z = 1 — 2 and then decreases again at 
higher redshifts. Using the results of a parallel submillimetre survey of the local 
universe, we show that both the submillimetre source density and background 
can be explained if the submillimetre luminosity density evolves in a similar 
way to the f/U-luminosity density. Thus, if these sources are ellipticals in the 
process of formation, they may be forming at relatively modest redshifts. 


^Department of Physics and Astronomy, Cardiff University, P.O. Box 913, Cardiff CF2 3YB, UK 

^Department of Astronomy, University of Toronto, 60 St. George Street, Toronto, Ontario MSS lAl, 
Canada 

^Mullard Space Science Laboratory, University of London, Holmbury St. Mary, Dorking, Surrey RH5 
6NT, UK 

"^CITA, University of Toronto, 60 St. George Street, Toronto, Ontario MSS lAl, Canada 
^DAEC, Observatoire de Paris-Meudon, 9219S Meudon, France 

®Laboratoire d’Astronomie Spatiale, Traverse de Siphon, B.P.8 13376, Marseilles Cedex 12, France 
^Dominion Astrophysical Observatory, S071 W. Saanich Road, Victoria, BC, V8X 4M6, Canada 



2 


1. Introduction 


The recent confirmation of a submillimetre background (|Puget et al. 199^ , [hixsen 
|et al. 19^ , |Hauser et al. 1995[ ) with approximately the same integrated intensity as the 
optical background ( pozzetti et ah 1998| , IVogeley 1998|) shows that to understand galaxy 
formation and evolution it is crucial to take into account that part of a galaxy’s starlight 
that is absorbed by dust and re-radiated at longer wavelengths. The integrated background 
is the repository of all the stellar energy that has ever been emitted (with a 1 + z weighting 
factor), and so the contribution of a particular class of galaxies to the background is a 
measure of its importance in the star-formation history of the universe. Determining the 
relative contribution to the background of different galaxy classes may be of particular 
importance for answering the question of how ellipticals and spiral bulges formed. The 
quantity of metals associated with nearby early-type galaxies (pdmunds and Phillipps 1997| ) 
implies that early-type galaxies should be responsible for about half the background (§4). 
Nevertheless, optical surveys have failed to hnd convincing evidence of proto-spheroids. 
Possible explanations of this failure are (a) that when the stars formed they were spread 
over a large volume of space before gradually coalescing to form galaxies (explaining the 
absence of discrete sources), (b) that the starbursts are hidden by dust. In both cases, 
these objects will still be contributing to the background; in the latter case they should be 
dominating the submillimetre background ([Bond et al. 1986] , [Dales and Edmunds 1997| ). 

The prospects for resolving the submillimetre background are good because of the 
commissioning of the SCUBA submillimetre array ([Holland et al. 1998] , [Gear et al. 1998[) on 
the James Clerk Maxwell Telescope. Last year a group carried out a submillimetre survey 
of two clusters, using the gravitational magnihcation by the clusters as an aid for Ending 
distant submillimetre sources. They found four sources and concluded that optical surveys 
had significantly underestimated the number of star-forming galaxies in the distant universe 
( [Smail et al. 1997| , SIB, [Blain et al. 199^) . We and others ( [Barger et al. 199^ , [Hughes ^ 


[al. 199(j[) have now started ‘blind-field’ surveys with SCUBA, which do not have the aid 
of gravitational amplification but do not suffer from the uncertainties of the cluster lens 
method. We are carrying out our survey in fields used for the Canada-France Redshift 
Survey ( [Lilly et al. 1995[) and our aim is eventually to survey 200 arcmin^ of the sky. In 
this paper we present the first results of a survey of an area of 22 arcmin^ and analyse these 
using the results of a parallel survey we are carrying out of the submillimetre properties of 
the local universe ( Punne et al. 1998[ ). A second paper ([Lilly et al. 1998] , hereafter Paper 
H) will discuss the optical counterparts to the distant submillimetre sources. 
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2. Observations and Data Reduction 


We observed three fields simultaneously at 450 and 850 /xm on 12 nights during 
January and March 1998, using the same observing procedure as SIB, with the exception 
that we chopped and nodded 30 arcsec EW. We divided our observations into integrations 
lasting about one hour, checking the telescope pointing before and after each integration. 
The optical depth during these nights varied between 0.57 and 1.65 at 450 /xm and between 
0.13 and 0.29 at 850 /xm. For the first field we simply made repeated observations at 
the same position (03^ 02™ 38.6'^ 00° 07' 40" [2000.0]), giving a total integration time of 
36.096 ks. For the second field we made eight integrations (21.504 ks in total) at one 
position (10^ 00™ 40.4® 25° 14' 20" [2000.0]) and four integrations (12.8 ks) at a second 
position (10^ 00™ 38.4® 25° 14' 20" [2000.0])) 27 arcsec away from the first. For the final 
field we made observations in a mosaic pattern such that each point in the final map 
was constructed from data received at nine locations within the instrument field of view. 
Since the instrument is fixed in Naysmith coordinates, and thus the bolometers move on 
the sky, this ensures that a point on the final map contains data from a large fraction of 
the available bolometers, reducing the possibility of problems with individual bolometers 
generating spurious sources. We reduced our data in the standard way using the SURF 
package ([Jenness 1997| , [Smail et al. 19971) . The final maps are shown in Figs 1 & 2. 


3. Source Extraction 

Since any real object is likely to be unresolved and will appear on a map as a positive 
source with two negative sources 30 arcsec on either side, we used template analysis 
( philips and Davies 1991]) to look for sources of this kind on the 850-/xm maps, using as 
a template a map of Uranus made with the same 30-arcsec chop. The negative emission 
is actually more complex than this because the chop direction is only precisely EW at the 
beginning of each integration, then remaining fixed at a constant azimuthal angle (this 
causes the slight rotation of the negative lobes visible in Fig. 2a), but this is a reasonable 
approximation. After trimming off the edge of each image, which is populated by noise 
artefacts due to incomplete sampling of the sky near the array’s edge, we made both a map 
of the correlation coefficient (CC), a measure of how closely the emission resembles the 
template, and a map of the convolution of the image with the template (Fig. 2). The noise 
on the template-convolved images is 0.56, 0.76, and 0.92 mJy per beam for the 3^, 10^, and 
14^ fields. We included a source in our catalogue only if both S/N > 3 on the convolved 
map and CC > 0.4, the latter criterion reducing the chance of noise spikes, which do not 
have the distinctive negative lobes, being included. 







4 


We checked the reality of our sources in two ways. First, we created the negative 
of each map and conhrmed that no sources were found using the same selection criteria. 
Second, we generated Monte-Carlo simulations of the helds, replacing the data for each 
bolometer (which is sampled every second) with the output from a gaussian random-number 
generator with the same standard deviation as for that bolometer (measured separately 
for each one-hour integration). The simulations showed that spurious sources tended to 
be concentrated towards the edges of our helds as expected, and we used the simulations 
to dehne ‘reliable regions’ for each map (Figs 1 & 2). The sources listed in Table 1 are 
from these regions, and the simulations suggest that, on average, one of these could be 
spurious. Because of the possibility that the rotation of the chop direction could affect 
huxes estimated by template-htting, we also estimated huxes by aperture photometry, 
hnding no signihcant difference for the two methods. 


Fig. 3 shows the integral source counts from our survey. An incompleteness analysis 
similar to that in SIB implies we are missing 0, 19, and 48% of the sources with huxes of 6, 
4, and 3 mJy respectively. Our source counts are consistent with those of the other small 
blank-held surveys ( [Barger et ah 199"^ , [Hughes et ah 199'8D and of SIB. 


Discussion 


The sources we have detected at 850 /rm account for 19% of the background emission 
at 850/im (Fig. 4), rising to 30% if a correction is made for incompleteness, and >10% 
of the background emission at 450/im. Our four measurements and seven limits on the 
450/rm/850/rm ratio are consistent with the average ratio being similar to that of the 
background emission, and thus consistent with the 850/im sources being representative 
of the background at shorter submillimetre wavelengths as well. On the assumption that 
the sub-mm objects detected in our survey make a negligible contribution to the optical 
background (Paper II), these objects account for ~11% of the total optical-near-IR plus 
submillimetre background ([Dwek et ah 1998f) , rising to ~17% if the incompleteness 
correction is made. 

It is possible that this level of background emission might be explained by AON. If 
the black-hole masses measured in nearby galaxies by HST ( [Ford et ah 1991^ ) are typical 
of galaxies in general, and if one makes the assumption that 10% of the mass accreting 
into a black hole is turned into energy, then the AGN background would be ~10% of 
the background from stars. However, given the uncertainties in this calculation, we will 
henceforth make the more conservative assumption that we are observing emission from 
dust heated by stars. 
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The integrated background emission produced by a population at a redshift z is related 
to the smoothed out cosmic density of processed material < p{Z + Ay) > produced by 
these objects by I ( [Fagel 19^) . Since our sources make up a 
signihcant fraction of the background, they must represent a population which has been 
responsible for a signihcant fraction of all the star formation that has ever occured. In Paper 
11 we show that these sources are similar to the galaxies with huge far-infrared luminosities 
discovered with IRAS ([Sanders and Mirabel 19961) . This combination of high individual 
bolometric luminosities and a signihcant contribution to the background emission is exactly 
what one expects for one mode of proto-spheroid formation (§1), and it is hard to resist the 
conclusion that these are these long-sought objects (see also Parger et ah 1998] , [Hughes ^ 
ah 1998[) . Note that, with their huge bolometric luminosities and implied star-formation 
rates, it has long been suspected that ultra-luminous IRAS galaxies might be ellipticals in 
the process of formation ( [Wright et al. 1990|) . What has been lacking is direct evidence 
that objects like this at high redshift are responsible for a signihcant fraction of the star 
formation in the Universe. 


An obvious question to ask is whether these objects are at z ~ 1 — 2, as are the galaxies 
that dominate the optical background (Pozzetti et ah 1998[) , or whether they are at higher 
redshifts, as one might suspect from the ages of the stellar populations in nearby ellipticals 
( Power et ah 1992 ). We will address this question by making the null hypothesis that the 
submillimetre luminosity density evolves in the same way as the UV luminosity density 
( [Madau et ah 19^ and predicting the submillimetre source counts and background. If 
the observed counts and background exceed our predictions, this might be evidence for 
a population of galaxies at z >> 1. In Paper 11 we will examine this question using the 
optical identideations. 

The biggest problem in modelling the submillimetre source counts and background 
has been the lack of information about the submillimetre properties of nearby galaxies. 
Without this information, one is forced to start from a local IRAS 60-/im luminosity 
function and make an assumption about the average far-infrared/sub millimetre SED of 
galaxies, an extrapolation of a factor of over 10 in wavelength (Pales and Edmunds 199^) . 
To overcome this problem, we have been carrying out a parallel submillimetre survey of 
the local universe (Punne et al. 19^) . By measuring submillimetre fluxes for galaxies 
drawn from statistically-complete samples of IRAS and optically-selected galaxies, we 
can both determine the SED’s for different types of galaxy and (using accessible-volume 
techniques— Avni and Bahcall 1980[) the local sub millimetre luminosity function. As 
yet our local 850-/im luminosity function is only well constrained at high submillimetre 
luminosities, which is acceptable for making source-count predictions, since it is only the 
high-luminosity sources which contribute at the flux densities of interest, but not for 































6 


estimating the background. Fig. 3 shows our source-count predictions. We implemented 
our null hypothesis in two ways: (i) by having the amplitude of the local luminosity function 
evolve with the same form as the luminosity density; (ii) by having the luminosities of the 
sources evolve with this form. Fig. 3 shows that with luminosity evolution, but not with 
space-density evolution, the predicted counts match the observed counts fairly well. If one 
assumes no evolution at all, the observed source counts exceed the predicted counts by a 
large factor, as noted by pmail et ah 1997 . 

For modelling the background we are thrown back on the IRAS 60-/im luminosity 
function ( [Lawrence et ah 1986| ), but we do have the advantage over previous studies of 
empirical results on the far-IR-submillimetre SED’s of galaxies. The 60, 100, and 850-/im 
fluxes of IRAS galaxies are well htted by a single-temperature dust spectrum with a dust- 
emissivity index of 1 and Td = 7.blogi^LpiR — 42.5 over the range 10.0 < IoqiqLfjr < 12.5, 
Lfir being dehned as in Lawrence et ah (Punne et al. 19^) . We do not claim that 
the dust-emissivity index is 1 and that there is dust of only one temperature in these 
galaxies, merely that this is a good empirical ht to the fluxes. Fig. 4 shows predictions 
using the IRAS luminosity function and this assumption about galaxy SED’s. The 
measured background is far in excess of the no-evolution model, but the model in which 
the submillimetre luminosity density evolves in the same way as the UV luminosity density 
matches the spectral shape and amplitude of the background remarkably well. 


Our conclusion that the null hypothesis of evolution like that in the optical waveband 
is acceptable (further supported by the analysis in Paper II) should be contrasted with 
recent claims to the contrary (plain et al. 199^ , [Hughes et al. 19^) . The difference is not 
in the data, but in the assumptions made about the local SED’s. Blain et al., for example, 
assume T^ = 38 K and n = 1.5. When compared with our assumptions about SED’s, this 
makes a factor of ~9 difference in intensity at 850 /im. 


The level of the background can also be used to argue against a major population 
of dust-enshrouded proto-ellipticals at z > 5. Let us assume that all the stars in this 
hypothetical population form instantaneously at a redshift zp and that the SED from one 
of these objects is the same as that of the archetypical local IRAS galaxy Arp 220 (the 
precise SED is not important). This is enough to dehne the shape of the background. The 
maximum amplitude of this background is then set by the observed CODE background (Fig. 
4). We then subtract this maximal elliptical background from the observed background to 
get a remnant background, which we assume is from the spiral/irregular population. If we 
then assume that late-type galaxies go through their period of maximum star formation at 
z ~ 1, we can then use the equation above and the ratio of the two backgrounds to predict 
the relative contributions to the present-day metal abundance of the two populations. If 
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zp = 5 the dusty ellipticals can produce <38% of the metals of the other population and 
if zp = 10 this falls to <8%. Observations suggest, however, that the ratio of the global 
metal abundance associated with ellipticals to that that associated with spirals/irregulars 
lies between 1 and 6 ( |Pagel 1997| , [Edmunds and Phillipps 1997|) . This discrepancy again 
suggests that ellipticals form at relatively modest redshifts. 


We are grateful to the many people that have contributed to the development of 
SCUBA and to the JAC staff for their excellent support during this program, in particular 
Wayne Holland for help in a multitude of areas. We also thank Rob Ivison for some useful 
advice. Loretta Dunne is supported by a PPARC studentship. Research by Dick Bond 
and Simon Lilly is supported by the Natural Sciences and Engineering Research Council of 
Canada and by the Canadian Institute of Advanced Research. The JCMT is operated by 
the Joint Astronomy Center on behalf of the UK Particle Physics and Astronomy Research 
Council, the Netherlands Organization for Scientihc Research and the Canadian National 
Research Council. 



Table 1. Source Catalogue. 


Name 

RA (2000.0) 

Dec (2000.0) 

S/N 

GG 

Ssso/im (mJy) 

8450^111 (mJy) 

CFRS03A 

03 02 36.27 

00 08 13 

6.1 

0.61 

3.4±0.6 

17±7 

CFRSlOA 

10 00 38.19 

25 14 50 

4.2 


4.8±1.0 

23±9 

CFRSlOB 

10 00 37.10 

25 14 58 

5.0 


4.7±1.0 

21±8 

CFRSIOC 

10 00 37.23 

25 15 12 

3.8 


2 . 8 ± 1.0 

<22 

CFRSIOD 

10 00 36.97 

25 14 42 

3.6 


3.9±1.0 

<22 

CFRS14A 

14 17 40.3 

52 29 08 

00 

bo 

0.78 

8 . 8 ± 1.1 

<31 

CFRS14B 

14 17 51.8 

52 30 29 

6.3 

0.64 

5.2±0.9 

<31 

CFRS14C 

14 17 33.6 

52 30 49 

3.4 

0.46 

3.8±1.1 

<31 

CFRS14D 

14 18 02.5 

52 30 53 

3.6 

0.47 

3.2±0.9 

<31 

CFRS14E 

14 18 02.6 

52 30 17 

3.6 

0.54 

3.5±1.0 

<31 

CFRS14F 

14 17 42.2 

52 30 27 

3.1 

0.43 

2.7±0.9 

20±11 

CFRS14G 

14 17 29.7 

52 31 15 

3.7 

0.45 

5.7±1.6 



Note. — Col. 1: name of source. Cols 2 & 3: position (2000.0). The positions were 
measured from the 850-/im maps, except for those of CFRS03A, CFRSlOA, and CFRSlOB, 
which were measured from the 450-/im maps. Col. 4: Signal-to-noise. This is the peak signal 
on the map produced by convolving the raw map with the template (see text) divided by 
the noise on that map, except in the cases of the 10 -hour sources, where the signal-to-noise 
was measured from the raw map. Col. 5: the correlation coefficient. The values for the 10^ 
sources are missing because template-fitting was not used for this field. Col. 6 : flux at 850^m. 
These were obtained by template-fitting, except those for the 10-hour sources were obtained 
by aperture photometry, using an aperture of diameter 20 arcsec. The errors include a 5% 
calibration error. Col. 7: fluxes at 450/im. These were obtained by aperture photometry, 
using an aperture of diameter 12 arcsec, except that the flux for CFRS03A was obtained by 
template-fitting. The errors include a calibration error of 20%. The limits are 3 ct upper limits. 
There is no measurement for CFRS14G because the slightly smaller field-of-view of SCUBA 
at 450 fxm means that at this wavelength it falls in a noisy edge region. 
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Fig. 1.— (Plate 1) Images at 850 /im of the 3^ (la) and the 10^ (lb) fields. The beam of 
the telescope has a FWHM of 14 arcsec at 850 /xm. The 3^ image has been convolved with 
a 14-arcsec ganssian (FWHM), the 10^ image with an 8-arcsec ganssian. The contonr lines 
are drawn at ±2, 3, 4, 5....a. The dashed line shows the region which is believed to be free of 
edge artefacts and to have reasonably nniform noise properties. This region is ~110 arcsec 
across for the 3^ field and ~120 across (in the vertical direction) for the 10^ held. 

Fig. 2.— (Plate 2) Images at 850 /xm of the 14^ held. The npper image (2a) has been 
convolved with a 14-arcsec ganssian, the lower image (2b) has been convolved with the 
sonrce template (see text). Both images have been divided by the sqnare root of a variance 
array prodnced from a Monte-Carlo simnlation of the held, which corrects for nneven noise 
across the images. The contonr lines are drawn at ±2, 3, 4, 5....a. As the convolntion 
with the template nses information from both the positive and negative ‘lobes’, there are 
more prominent sonrces in 2b than in 2a, althongh one has to beware spnrions sonrces being 
created 60 arcsec EW from a bright sonrce (e.g. sonrce Z). The dashed line in the lower 
image shows the region which is believed to be free of edge artefacts. This is ~390 arcsec 
across in the horizontal direction. 


Fig. 3.— Integral sonrce connts at 850/im from onr snrvey (dots) and other surveys: triangle 
( [Hughes et al. 1998|) , circle ( [Barger et al. 19'9B[ ) and cross ([Small et al. 1997| ). The 
errors on our counts are Poisson errors and are not independent. No correction has been to 
any of the counts for incompleteness, although the incompleteness for the other two blind- 
field surveys is claimed to be low. The predictions were made as described in the test. 
The key is as follows: no evolution, Hq = 1—squares; density evolution, Hq = 0—dots in 
circles; density evolution, Hq = 1—circles; luminosity evolution, Hq = 0—stars; luminosity 
evolution, Hq = 1—triangles. 


Fig. 4.— The dashed line and the points show, respectively, the submillimetre background 
estimated from the results of the FIRAS ([Fixsen et al. 199'8[) and DIRBE (pdauser et al. 1998[) 
experiments on COBE. The two lower limits show the limits on the background produced 
by simply adding up the number of sources that we observe. The solid line and the dotted 
lines show backgrounds predicted using an IRAS luminosity function and information on 
the SED’s of galaxies from our nearby galaxy survey. The solid line is for no evolution, the 
dotted line for a model in which the submillimetre luminosity density evolves like the UV 
luminosity density. The two other lines are the backgrounds predicted for a model in which 
all the stars in ellipticals form instantaneously at a redshift zp^ with the stars being hidden 
by dust. The predicted backgrounds have been normalized so that they are the maximum 
backgrounds consistent with the FIRAS measurement. The higher curve is for zp = 5, the 
lower for zp = 10. 
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